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ABSTRACT 



We derive a simple analytical solution for the evolution of a close binary 
with nuclear time-scale driven mass transfer from a giant. This solution is 
based on the well-known fact that the luminosity and the radius of a giant 
scale to a good approximation as simple power laws of the mass Mc of the 
■ degenerate helium core. Comparison with results of numerical calculations 

I by Webbink, Rappaport & Savonije (1983) show the analytical solution and 

the power law approximation to be quite accurate. The analytical solution 
presented does also allow (in parametrized form) for non-conservative mass 
. transfer. Furthermore it is shown that the near constancy of the mass transfer 

' rate over most of the mass transfer phase seen in the results by Webbink, 

0^ ■ Rappaport & Savonije is not a generic feature of this type of evolution but 

rather a consequence of a particular choice of parameters. The analytical 
(~| . solution also demonstrates that the level of mass transfer is largely set by the 

Oh' core mass of the giant at the onset of mass transfer. Finally we show that the 

Q . model is selfconsistent and discuss its applicability to low-mass X-ray binaries. 

Key words: binaries: close - stars: evolution - stars: pulsars: general - 
I X-rays: stars 

X ■ 1 INTRODUCTION 

Stable mass transfer on to a (compact) accretor which is driven by the nuclear evolution 
of a giant is a process of considerable astrophysical interest. This type of evolution was 
first encountered in numerical simulations by (among others) Kippenhahn, Kohl & Weigert 
(1967) and then studied in some detail by Refsdal & Weigert (1969, 1971) in the context of 
the formation of low-mass white dwarfs in binaries. Later Webbink, Rappaport & Savonije 
(1983, hereafter WRS) and, independently, Taam (1983) have proposed that some of the 
most luminous among the galactic low-mass X-ray binaries could be powered by nuclear 
time-scale driven mass transfer from a giant. Finally with the detection of the millisecond 
pulsar binary PSR 1953-1-29 by Boriakoff, Buccheri & Fauci (1983) it became immediately 
clear that this system (and many other similar ones which have been found since then) must 
be the end product of an evolution of the type described by WRS and Taam (1983), see e.g. 
Joss & Rappaport (1983), Paczynski (1983), Savonije (1983), and Rappaport et al. (1995, 
hereafter RPJDH) for a more recent reference. 

Whereas Taam (1983) performed full stellar structure calculations WRS took a simpler 
approach by making use of the well-known fact that giants obey a core mass-luminosity 
and a core mass-radius relation. WRS derived these relations from full stellar equilibrium 
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models of giants and approximated them as power series of the mass of the degenerate hehum 
core Mc- By using the core mass-luminosity and core mass-radius relation calculating the 
evolution of a binary with nuclear time-scale driven mass transfer from a giant reduces to 
preforming a simple integral in time which, however, in the approach taken by WRS has still 
to be performed numerically. The purpose of this paper is first to show that the evolution 
of such a binary can be described by an analytical solution if the core mass-luminosity and 
core mass-radius relation are simple power laws of M^, and second to show that the results 
of full stellar structure calculations of giants given by WRS can indeed be approximated 
with little loss of accuracy by simple power laws. Moreover, the analytical solutions which 
we present are more general than the numerical ones of WRS in the sense that we allow for 
non-conservative mass transfer and a non-zero mass radius exponent of the donor star (see 
below section |^). 

The advantages of having an analytical solution as compared to a numerical one are 
immediatedly clear: First, if sufficiently simple (as in our case), such solutions provide us 
with direct physical insight into how this type of evolution works, i.e. which parameters 
determine the characteristc properties and which ones are of minor importance. Second, 
whereas a numerical solution has to be computed separately for each set of input parameters 
of interest, an analytical solution provides us with the whole manifold of solutions, i.e. we get 
the explicit dependence of the solution on the input parameters. Among other things, this 
allows us to examine whether certain properties found in numerical solutions are generic or 
a simple consequence of a particular choice of parameters. A case in point is the well-known 
property of the solutions presented by WRS, namely that the mass transfer rate stays almost 
constant over most of the mass transfer phase. 

Another area of application in which simple analytical solutions are superior to numerical 
ones is population synthesis. This is because population synthesis calculations usually involve 
integrals over the complete manifold of solutions of the particular evolutionary phase under 
consideration. Using numerical solutions for that task requires computation of a sufficiently 
dense grid of evolutions over a multi-dimensional parameter space and such computations 
can easily become prohibitively expensive (see e.g. Kolb 1993 for a detailed discussion in the 
context of cataclysmic variables). Using analytical solutions instead is comparatively easy 
and cheap. 

This paper is organized as follows: In section |] we present the basic model assumptions 
in the framework of which we shall work, in section |] we give an analytical solution for 
the nuclear evolution of a giant, and in section | we derive expressions for the mass radius 
exponent of the donor star's critical Roche radius. Finally, in section |] we present the 
analytical solution for nuclear time-scale driven mass transfer from a giant. In order to 
compare our results with those of WRS we derive in section ^ the parameters of the power 
law approximations of the core mass-luminosity and core mass-radius relations using the 
data from WRS and RP JDH. In section ^ we then compare the results obtained from the 
analytical solutions with the corresponding numerical ones of WRS. In section ^ we examine 
whether the model is selfconsistent and whether it is applicable to low-mass X-ray binaries. 
A discussion based on the analytical solution of the key features of nuclear time-scale driven 
mass transfer from a giant follows in section |^, and our main conclusions are summarized in 
the final section PH. 



2 BASIC MODEL ASSUMPTIONS 

In the following we shall list the basic model assumptions in the framework of which we 
shall work. These assumptions are: 
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(i) The donor star (of mass M2) is on the first giant branch and has a degenerate hehum 
core of mass M^.. Nuclear luminosity comes from hydrogen shell burning exclusively. 

(ii) M2 is small enough that the binary is adiabatically and thermally stable against mass 
transfer. For this to hold we must have (e.g. Ritter 1988) 

Cad - Cr,2 > (1) 

and 

Ce - Cr,2 > 0, (2) 

where 

^^'={^^)^, 

is the adiabatic mass radius exponent of the donor star, 
. fd\nR2\ 

the corresponding thermal equilibrium mass radius exponent, and 

the mass radius exponent of its critical Roche radius (see e.g. Ritter 1996 for details). In 
(I) the subscript * indicates that for computing this quantity one has to specify how, as 
a consequence of mass transfer, mass and orbital angular momentum are lost from and 
redistributed within the binary system (see below section^). 

For stars on the first giant branch and with small relative core mass M^/M Cad ~ —1/3. 
Cad increases with relative core mass and becomes positive if M^/M ^ 0.2 (Hjellming & 
Webbink 1987). Therefore, as we shall see, adiabatic stability of mass transfer is less critical 
than thermal stability. 

Assuming that the donor obeys a core mass-luminosity relation, i.e. that L = L{Mc) and 
OL/dM = 0, where L is the (nuclear) luminosity, Ce follows from the Stefan-Boltzmann law: 

^ 1 ( d\nL \ _ ( d\nT,s \ ^ ( d\nT,A 
2\dlnM)^ [dlnMj^ [dlnMJ^ 

It is well known that a giant, as long as it has a sufficiently deep outer convective envelope, 
in the HRD stays close to and evolves along the Hayashi line to higher luminosity. This 
property allows us to get an approximation for Cc- Approximating the Hayashi line in the 
HRD in the form (e.g. Kippenhahn & Weigert 1990) 

logL = alogTcff + 61ogM + c (7) 
this together with (H) yields 

Ce ~ -2b/ a . (8) 

Numerical calculations typically yield —0.3 ;^ Cc ~ —0.2. Thus, except for small relative core 
mass, Ce < Cad, i-G. the criterion for thermal stability (||) against mass transfer is stronger 
than the one for adiabatic stability (|I]). 

(iii) Mass transfer is driven only by the donor's nuclear evolution, i.e. we assume that 
there is no systemic angular momentum loss in the absence of mass transfer. Consequential 
angular momentum loss (see e.g. King & Kolb 1995 for a discussion) is, however, not excluded 
a priori. 
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(iv) Mass transfer is so slow that the donor star remains close to thermal equilibrium. 
Accordingly we assume that its radius R2 and luminosity L2 are given by the corresponding 
thermal equlilibrium values, i.e. R2 = R2,e and L2 = L2^e, and hence Cefr = {dlnR2/dlnM2) = 
Ce- This condition is well fulfilled because mass transfer occurs on the nuclear time scale Tnuc 
which is much longer than the Kelvin-Helmholtz time tkh = GM2^/-R2-^2- We shall later 
(section ^ verify the validity of this assumption. 

With the above assumptions we can now write the mass transfer rate, i.e. the donor's 
mass loss rate as (e.g. Ritter 1996) 



Ce - Cr,2 V 9t 

where (91n/22/<9t)nuc is the inverse time-scale of the donor's expansion due to nuclear evo- 
lution. Next we shall work out in more detail ((91ni?2/<9t)nuc and Cr,2- 



3 NUCLEAR EVOLUTION OF A GIANT 

Here we are going to exploit the fact that to a very good approximation the luminosity L of 
a star on the first giant branch scales as a simple power law of the mass of the degenerate 
helium core but does not depend on its total mass (Refsdal & Weigert 1970). Accordingly 
we make the following ansatz for L: 

L(M.)=Lof^V (10) 



M, 



0- 



It then follows from Eqs. (^ and (|^ that the radius of such a star cannot at the same time 
depend only on the core mass Mc but must also depend on its total mass M. Because along 
the Hayashi line the effective temperature remains almost constant it follows from (ITO) and 



the Stefan-Boltzmann law that the appropriate form of the 'core mass-radius relation' is 

Here, Lq, A, Rq and p are parameters which have to be determined from full numerical calcu- 
lations. We shall give examples (taken form WRS and RP JDH) in section ^. Using homology 
techniques Refsdal & Weigert (1970) could derive an expression for A, the numerical value 
of which depends on the conditions in the hydrogen shell source. A typical value for Pop. I 
giants is A ~ 7 — 8. 

Hydrogen burning in the shell source around the degenerate core adds to the core mass 
at a rate Mc. Thereby the nuclear luminosity generated is 

L = Lnuc = XQM, , (12) 

where X is the hydrogen mass fraction in the star's envelope and Q ~ 6 10^^ erg g^^ the net 
energy yield of hydrogen burning per unit mass. In writing ([T2|) we have assumed that the 
star is in thermal equilibrium, i.e. that L = Lnuc or, for practical purposes, that \Lg\ <^ -^nuc, 
where Lg is the gravo-thermal luminosity. 



(jTO) and (|T2|) together yield a simple differential equation for Mc{t) with the solution 



Mc(t) = Me(t = 0) ( 1 - ) ' , (13) 



where 
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(A - l)Lo 



"Mc(t = 0) 

Mo 



l-A 



(14) 



is the time over which, formally, the core mass reaches infinity. It is also the characteristic 
time for nuclear evolution on the first giant branch. (13) combined with respectively (10) or 
(0) yields 



m 

and 

m 



M,{t = 0) 



M, 



© 



1 - 



A 
l-A 





1' 


"Mc(t = 0)' 








[ Mo J 


too) 



(15) 



(16) 



Here Mc{t = 0) is the core mass at the beginning of the evolution on the first giant branch. 
This quantity has also to be taken from full stellar evolution calculations or to be treated 
clS db free parameter. 

From (|T6|) we finally obtain 

n 1 

(17) 



'dlnR' 



Mr 



P 



dt 



A-1 t 



t 



the expression to be inserted in 



4 COMPUTING Cr,2 

Here we follow standard arguments which have been given many times in the literature (e.g. 
Soberman, Phinney & van den Heuvel 1997 and references therein). Therefore, we shall give 
here only the definitions we use and the final results. In the context of this paper, we assume 
that mass transfer occurs from the secondary and that the primary (of mass Mi) accretes 
(on average) a fraction rj of the transferred mass, i.e. 

dMi = -r]dM2 (18) 

and that the mass d{Mi + dM2) = (1 — r])dM2 leaves the binary system carrying away the 
orbital angular momentum 

d(M, + M.) 

M1 + M2 ' ^ ^ 

where, for the moment at least, rj and u are free parameters, 

is the orbital angular momentum, and a the orbital separation. Writing for the secondary's 
critical Roche radius 

i?R,2 = a/2(g), (21) 
where 

g = ^ (22) 
is the mass ratio, one derives 
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Ca. = (l-.,)^ + ^-2-Afl + ^), (23) 

1 + g q \ l) 

with 

A = (24) 

dmq 

As a further simphfication we use for /2 in (^) and (|2^ ) the approximation (Paczynski 1971) 

^-^-^^ ^^^^ 

which yields 

/32(g) = . (26) 

^ 3(1 + g) ^ ' 

Inserting (^6]) in (|23D finally yields 

(l-r/)(2z/+l) 2Z/ , r^ + g . . 

1 + g g 3(1 + g) 

In the following we shall also give the results for two specific choices of the parameters r] 
and V. 

4.1 Conservative mass transfer 

Here r] = 1 and z/ = 0. From (E3) and (PB|) we then have 



Cb. = ^. (28) 

3g 

4.2 Isotropic stellar wind from the primary 

Here we are considering the case that all the transferred mass is ejected from the primary 
in an isotropic stellar wind carrying with it the specific orbital angular momentum of the 
primary. Accordingly 77 = and v = 1/q. Inserted in (^) together with (|26| ) this yields 

6 - 5g2 - 3g 

'''' 



5 EVOLUTION WITH MASS TRANSFER FROM A GIANT 

Now that we have determined (9 In i?/9t) nuc and Cr,2 in equation (^ we can solve that 
equation. First we are going to give a rather general solution before discussing special cases 
of interest. 

5.1 General analytical solution 

General solution in this context means that we allow for a constant value rj > 0, either a 
constant value of > or z/ = 1/g, and a constant value of Ce subject to the condition @. 
Inserting ([T7|) and ( ^31) with (|26|) in @ yields an ordinary differential equation for M2{t), 
or rather, as we shall see, for t{M2) which is solved by separating the variables M2 and t: 

Jm,, m2 ^0 V dr J 
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The solution of (130|) can be written as follows: 



t{M2) = tc 



'Mr 



-P2 



-P3 



-Pi ■ 



(31) 



where the exponents pi, P2, P3 and p4 are given in Table 1 for various cases of interest. With 



(0) the mass transfer rate can be written as 

A-l 

I n/,n //!//„; \ 

M2 



Mo 



t 



1 



Ce-CR,2 VM©; vmq, 

Inserting (^1]) in ([T5[) and (|l^) finally also yields 



Lo — Lf] 



and 



Ro = Ro 



' Me/ 



A r 



'M^ 



pi 



M2Y' 



M2 



Cc 



M2A. 



Ml 



' Ml + M2 ^ 
Mi,i + M2i 



P3 



jv:( 2 i ~ ^^-^ 2 



Pi 



Mo 



P2 



Ml + M2 



P3 



2 i ~ 2 



(32) 



(33) 



(34) 



Here Mi ;, M2,i and M^^i are respectively the initial values, i.e. taken at t = 0, of Mi, M2, 
and Mc- Note that from ([T8|) we have 



Ml = Mi,i + 7] (M2,i - M2) . 
As can be seen the equations 



(35) 

I) (together with (|3T|) ), ( p3D and (|3^ do not depend explicitly 

on time t. The solutions are completely fixed by the parameters of the core mass-luminosity 
relation (|T0|) , i.e. Lq and A, the 'core mass-radius relation' (|ll|), i-e. -Ro, P and Ce, the initial 
values and the current values of both masses and the secondary's core mass. Kowing M2, 
i?2 and L2 as explicit functions of the initial values and current values of Mi, M2 and Mc 
other quantities of interest such as the orbital period P, the orbital separation a, the size 
of the accretor's Roche radius Rr^i and with it the size of the disc around the accretor 
and the secondary's effective temperature Teff,2 follow directly. What this solution does not 
provide explicitly is the total duration of the mass transfer phase., i.e. the time tf when 
the secondary's envelope is exhausted or, which is equivalent, the final core mass Mc^f. A 
very good estimate of the latter together with tf can be obtained by inserting ( PT| ) with 



M2 = M2,f ^ M^^i in 



Note that the resulting equation for M^^f cannot in general be 



solved in closed form. Reinserting the solution of this equation in (^ip finally yields also tf. 



5.2 Special cases 

Table 1 lists the powers Pi, P2, Ps and p^ which appear in equations ( pT]) - p^ ) for a number 
of special cases of interest. 

Case 1 is the most general case for constant values of rj, v and Ce- In case 2 the special 
choice V = 1/q with 77 > is made, whereas case 3 corresponds to case 2 in the limit = 0. 
Cases 4 and 5 correspond to conservative mass transfer with respectively Q = const. 7^ 
and Ce = 0. Case 6 finally corresponds to case 3 with = 0. 

Case 3 (or 6) corresponds to the situation in which the accretor ejects all the transferred 
mass with the specific orbital angular momentum of the accretor's orbit. This approximates 
what apparently must happen in systems in which the mass transfer rate to a neutron star 
exceeds the Eddington accretion rate by many orders of magnitude (see e.g. King & Ritter 
1999 for a discussion). 

We shall later mostly make use of case 5 because this is the approximation within which 
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Table 1. Parameters and coefficients entering the analytical solution 



case 


V 


V 


Cc 


Pi 


P2 




Vi 


P4 




1 


> 0, const. 


const. 


const. 


2(A-1) 


(3Cc+5){A- 


1) 


2(3i'+l){A-l) 







P 


3p 




3p 




2 


> 0, const. 




const. 


2(A-1) 


(3Cc+5){A- 


1) 


4{A-1) 









3p 




3p 




3 





i/g 


const. 





(3Cc+5){A- 
3p 


1) 


4{A-1) 
3p 


2(A- 
P 


1) 


4 


1 




const. 


2(A-1) 


(3Cc+5){A- 


1) 












P 


3p 






5 


1 







2(A-1) 
P 


5(A-1) 
3p 












6 













5(A-1) 
3p 




4{A-1) 
3p 


2(A- 
P 


1) 



WRS have performed their numerical calculations and with which we wish to compare the 
results of our analytical solution. Before we can do that we have to provide the parameters 
describing the core mass- luminosity and core mass-radius relation. 



6 APPROXIMATIONS FOR THE CORE MASS-LUMINOSITY AND 
CORE MASS RADIUS RELATION 

In general the values of Lq, A, Rq, p and Ce in equations ( ]TI1| ) and (|TT]) have to be determined 
from full numerical calculations of stars in thermal equilibrium on the first giant branch. 
Here we restrict ourselves to deriving approximate values of these parameters by using 
results given by WRS (their equations 4-6). Because WRS have assumed Ce = their 
relation (5) is a true core mass-radius relation in the sense that the radius of the star does 
not depend on the total mass M. Their results for Pop. I and Pop. II giants are shown 
graphically respectively in Fig. 1 and Fig. 2. It is immediately apparent that in the range 
O.15M0 ^ Mc ^ O.4M0 the relations L[M^ and R{Mc) are very nearly simple power laws 
in Mc- Based on the results of Refsdal & Weigert (1970) this is of course not surprising. 
Because our analytical solutions of section ^ are based on the power law approximations 
(p!0|) and ( p^D we fit the results of WRS shown in Figs. 1 and 2 accordingly. The resulting 
parameters, together with the chemical composition of the underlying numerical models 
(characterized in the usual way by the hydrogen mass fraction X and the the metallicity Z) 
are summarized in Table 2. The corresponding relations are shown as thin full (L) or dashed 
{R) lines in Figs. 1 and 2. 

For comparison we provide in Table 2 also approximate parameters of a power law fit in 
the range O.2M0 ^ ^ O.SMq to the core mass-radius relation derived by RPJDH (their 
equation 5). 



7 NUMERICAL RESULTS VERSUS THE ANALYTICAL MODEL 

Here we wish to compare the results obtained from the analytical approximation, i.e. from 
equations ( pij ) - (|5^), with the numerical results obtained by WRS. For this we have to insert 
the parameters listed in Table 2 into the expressions for pi, p2, ps, and p^ corresponding 
to case 5 in Table 1 and in equations (^TJ) - (|5). Rather than listing all the resulting 
equations for both Pop. I and Pop. II composition, we restrict ourselves to giving only the 
corresponding expressions for the mass transfer rate. With Q = 6 10^^ erg g~^ this yields 
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log Ten (K) 



3.60 



3.50 



Me 

Figure 1. Core mass-luminosity relation (full line), core 
mass-radius relation (long-dashed line) and core mass- 
effective temperature relation (short-dashed line) for Pop. I 
giants according to numerical calculations by WRS. Our sim- 
ple power law fits are shown respectively as thin full and thin 
long-dashed lines. The corresponding fit parameters are listed 
in Table 2. 




log Te„ (K) 



log 



Mc 



Figure 2. The same as Fig. 1, but for Pop. II giants. 



Table 2. Approximate parameters of the core mass-luminosity and core mass-radius relations derived from numerical results 



Pop. 


X 


z 


logi| A 




p 




source 


I 


0.70 


0.02 


6.3 8 


4.1 


5 





WRS 


II 


0.80 


0.0001 


4.5 5 


2.7 


3 





WRS 


I 


0.70 


0.02 




3.58 


4.21 





RPJDH 


II 


0.70 


0.001 




3.36 


4.17 





RPJDH 


-M2 


= 2.64 10" 






V 




Ml 


^O.25M0 


J 5Mi 


- 6M2 



Mo 



for Pop. I, and 
- M2 = 1.40 lO^^Moyr"^ 



^0, 



Ml 



Mo 



M2A, 



Ml,, 



MA'""" (MA-"" 



O.25M0/ 5M1-6M2 \Mq \M2i 



Ml, 



(36) 



(37) 



for Pop. II chemical composition. 
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I ' ' ' ' ' ' 1 

1 2 3 4 5 6 

t(10^ yr) 

Figure 3. Comparison of the result obtained with the ana- 
lytical solution (full line) with that of a matching numerical 
one by WRS (dashed line). Shown is the mass loss rate from 
a Pop. I giant as a function of the time for the parameters 
used by WRS in their figure 3, i.e. X = 0.70, Z = 0.02, 
Mi,i = 1.4Mq, M2,i = IMq, and Md = 0.26Mq. 



I ^ ^ ^ ^ ^ ^ 1 

1 2 3 4 5 6 

t(10^ yr) 

Figure 4. Comparison of the result obtained with the ana- 
lytical solution (full line) with that of a matching numerical 
one by WRS (dashed line). Shown is the mass loss rate from 
a Pop. II giant as a function of the time for the parameters 
used by WRS in their figure 4, i.e. X = 0.80, Z = lO"", 
A/i^i = 1.4Mq, M2,i = O.9M0, and M^-^ = CSIMq. 



Results of detailed numerical computations by WRS are shown in their figs. 3 and 4. 
Comparing the results obtained from our analytical solution for the same parameters with 
those of WRS shows that the differences are small, too small in fact to be adequately shown 
on the scale of these figures. That the differences in L{t) and R{t) are small was to be 
expected based on the quality of the fits used (see Figs. 1 and 2). That the differences in 
— M2 are small too is shown in Figs. 3 and 4. In that context we should emphasize that our 
fits of the L{Mc) and R{Mc) relations of WRS with power laws are the simplest possible 
ones (with integer powers A and p) and that no attempt has been made to adjust the power 
law parameters such as to minimize the differences between the numerical results of WRS 
and the analytical solution. With —M2 being in close agreement with the numerical results 
it follows that also M2{t) and P{t) should agree very well, which is indeed the case. Finally 
also Mc{t) must be in close agreement with the numerical result of WRS, otherwise the 
temporal evolution of the other quantities would not match either. 

In fact we can summarize the result of this comparison as follows: if the power law fits to 
the core mass-luminosity and core mass-radius relations obtained from full stellar structure 
calculations are adequate, the analytical model provides an adequate approximation to the 
results of numerical computations. 



8 SELFCONSISTENCY AND APPLICABILITY OF THE MODEL 
8.1 Deviation from thermal equilibrium 

The equations we have derived in section |^ are valid only to the extent that the donor star 
is not significantly driven out of thermal equilibrium as a consequence of mass loss. That is 
to say that the time-scale of mass loss 

dt 

must be significantly longer than the shortest time-scale on which the donor's convective 
envelope can react, i.e. longer than the Kelvin-Helmholtz time of the convective envelope 
(e.g. King et al. 1997a) 



© 0000 RAS, MNRAS 000, 000-000 



Analytical solution for the evolution of a binary with stable mass transfer from a giant 11 



3 Mee 

TKH,ce — TkH — ^ 



3 M - Mc GM^ 



M RL - f^'" 

Here M^e = M — is the mass of the convective envelope and tkh the Kelvin-Helmholtz 
time of the whole star. A measure for the star's deviation from thermal equilibrium is thus 
the the ratio TMa/TKH.ce- Using equations and ( |T5| ) 
write 

l-A 



r?!) in ( |38D and (^) we can now 



''"M2 — (Ce — Cr,2 

and 



Tkh,, 



3 GM©^ Ma - Mc / Ma 



2-Co 



7 -Ro-^o 
and, therefore, 



M, 



TM2 

TkH,( 



(Ce - C 



R,2j 



7 XQi?o M2 



to 



Mo. 



Mo 



-(A+P) 



1 

tr^ 



A+p 
A-1 



Cc-2 



3 pGMq M2 -Mc VM, 



0. 



M, 







1 - 



tr. 



(40) 



(41) 



(42) 



Inserting now for illustrative puposes the parameters of Table 2 which correspond to the 
results of WRS into (H) yields 



= 3.16 10^ 

TKH,cc 

for Pop. I, and 
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0. 



0.25 Mo 



0.25 Mo 



1 - 



t 



1 



(43) 



(44) 



for Pop. II chemical composition. 

What (1^) - (1^) demonstrate is that for all practical purposes TMj/TKH.ce ^ 1 unless the 
mass ratio is such that the binary system is close to thermal instability, i.e. that Ce~CR,2 "C 1. 
Thus mass transfer driven by nuclear evolution of a star on the first giant branch is always 
sufficiently slow for the donor star to remain essentially in thermal equilibrium. To this 
extent, the analytical model presented in section ^ is self consistent. 



8.2 Irradiation of the donor star 

Since one of the main areas of application of this model is the evolution of luminous X-ray 
binaries (see e.g. Taam 1983 and WRS) and their transformation into millisecond pulsar 
binaries (see e.g. RPJDH and references therein) we have to examine whether this model 
remains applicable if irradiation of the cool donor star by X-rays from the vicinity of a 
very compact accretor is taken into account. As has first been shown by Podsiadlowski 
(1991) cool stars, when exposed to sufficiently strong external heating by X-rays, respond 
by expansion on the thermal time-scale of the convective envelope, thereby giving rise to 
very high mass transfer rates in a semi-detached binary. In the context of mass transfer 
from a giant to a compact accretor the stability of mass transfer with taking into account 
irradiation has been examined in some detail by King et al. (1997a). Their finding was that 
if the accretion rate on to the compact star equals the mass transfer rate at any time, the 
evolutionary model discussed here is not in general applicable, because mass transfer is very 
likely to be unstable. Rather than being continuous, mass transfer occurs in cycles which are 
characterized by short high states during which very high mass transfer rates are achieved 
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and long low states during which the binary is detached. As has been shown by King et al. 
(1997a) this instability occurs mainly because for giants the ratio of mass loss time-scale 
to Kelvin-Helmholtz time (equation |2) is so large. It is therefore somewhat ironic that 
self-consistency of the model of nuclear time-scale driven stationary mass transfer from a 
giant is possible only at the price that such systems are in danger of being unstable to 
irradiation-driven mass transfer. 

However, the above-mentioned condition, namely that the accretion rate equals the mass 
transfer rate at any time, is not likely to be fulfilled. Rather, as has been shown by King 
et al. (1997b), the accretion disc in such binary systems is likely to be unstable and that 
even for the rather optimistic assumptions made by King et al. (1997b) about the efficiency 
of X-ray irradiation of the accretion disc. Based on newer results of Dubus et al. (1998) 
one could conclude that the accretion discs in essentially all systems with a giant donor are 
unstable. If this is the case, then accretion on to the compact star is spasmodic because 
the disc undergoes a dwarf nova-like limit-cycle instability. Estimates of the properties of 
the outbursts (e.g. King & Ritter 1998) show that they are characterized by a very small 
duty cycle. This, in turn, is now of great importance for the stability of mass transfer from 
the giant in the presence of X-ray irradiation: King et al. (1997a) show that if accretion on 
to a neutron star or a black hole is spasmodic with a small duty cycle, then mass transfer 
from the giant is stable despite the donor being irradiated. If that is the case, the model 
considered by Taam (1993), WRS and us is applicable. 

On the other hand, if the accretor is a white dwarf the situation is more complicated: If 
the core mass of the donor is ^ 0.2Mq, irradiation cannot destabilize mass transfer and, 
therefore, the model is applicable. Yet if ^ O.2M0, mass transfer is unstable. Although 
spasmodic accretion has a stabilizing effect, for mass transfer to become stable in this case 
the required small duty cycle has to be smaller still by a factor of ~ 10 ■ ■ ■ 100 than in the 
case of a black hole or a neutron star (see King et al. 1997a, fig. 4). Because the value of 
the duty cycle does not depend explicitely on the nature of the accretor but rather on the 
average mass transfer rate and on the size of the accretion disc, it appears unlikely that the 
duty cycle in a white dwarf system is systematically smaller than in a black hole or neutron 
star system. Therefore it is more likely that in a white dwarf system in which M^. ^ 0.2Mq 
mass transfer is unstable despite spasmodic accretion. In this case the model in question is 
not applicable. However, we should emphasize that detailed numerical calculations of mass 
transfer from a giant with taking into account the effects of irradiation of the donor star and 
of disc instabilites have yet to be done. Before results of such computations are available, a 
final decision about whether this model is applicable cannot be made. 



9 DISCUSSION 

Having established that the analytical approximation is both accurate and self consistent, 
and that the model is probably also applicable to the most important case of interest, i.e. 
to low-mass X-ray binaries, we can now proceed to discuss the properties of these solutions. 
One of the remarkable properties of this type of evolution which has already been pointed 
out by WRS is that the mass transfer rate remains almost constant for most of the mass 
transfer phase. Here we shall examine whether this is a generic property of mass transfer 
from a giant or whether this is due to a particular choice of parameters. Other aspects which 
we wish examine are the dependence of the temporal evolution on the initial parameters, in 
particular on the initial core mass Mc^i, and systematic differences between Pop. I and Pop. 
II systems. 
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Figure 5. Mass transfer rate computed from the analytical solution for three different cases: The full and dashed lines show 
respectively the mass transfer rate as a function of time in an evolution of a Pop. I and a Pop. II system with same initial 
parameters as those used in Figs. 3 and 4, but now for non-conservative mass transfer characterized by = and u = 1/q. The 
dash-dotted line shows the mass transfer rate in a conservative evolution of a Pop. I system with the parameters Mi j = SMq, 
M2,i = 2M0, and Mc,i = O.26M0. 



9.1 The mass transfer rate as a function of time 

One of the remarkable properties of the evolutions shown by WRS in their figs. 3 and 4 (see 
also our Figs. 3 and 4) is that the mass transfer rate remains constant within a small factor 
2) in the case of conservative mass transfer, and, as Fig. 5 shows, within an even smaller 
factor 1.5) in the case that rj = and u = 1/q throughout the evolution. This raises the 
question why this is so and whether this is a generic property of this type of evolution. For 
our discussion we combine and ( P!7D to write the mass transfer rate as 

_M,._i4 ^. (45) 

Cc — Cr,2 a — 1 too — t 



From ( [45D it is immediately seen, that it contains two counteracting factors: first, with time 
Ce — Cr,2 increases because q increases. This factor has thus the tendency to reduce —M2 
with time. On the other hand, t^o — t decreases with time and thus has the tendency to 
increase — M2. Therefore, there is a possibility that —M2 goes through a minimum during 
an evolution. Yet, an evolution does not necessarily include this minimum. Rather, whether 
or not —M2{t) has a minimum depends on the initial mass ratio qi and on the subsequent 
change of q with time. Starting from equations ( ^I]) and ( ^2]) and working out M2 = leads 
to a (rather complicated) equation for the mass ratio go at which —M2 becomes extremal. 
From the above discussion we know already that this extremum is a minimum. Denoting 
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the final mass ratio by gf = Mif/M2^f ~ Mi^f/M^^f, we can now distinguish three different 
cases: 1) if gcrit < Qi < % < Qf then the mass transfer rate goes through a minimum, 
2) if gcrit < Qo < Qi < Qi, the mass transfer rate will increase monotonically, and if 3) 
?crit < < Qf < Qo, the mass transfer rate will decrease monotonically. Here gcrit is the 
critical mass ratio for both thermal and adiabatic stability against mass transfer. Its value is 
given by the larger of the positive solutions of the equations Cad — Cr,2 = and Ce — Cr,2 = 
(see e.g. Ritter 1988). 

For illustrative purposes we give the values of gcrit and go for the cases we show in Figs. 
3-5. For these cases we have assumed with WRS Ce = and implicitly Ce < Cad- From this 
and (^) it follows that gcrit = 1.2 in the case of conservative mass transfer, and with (^Of ) 
that gcrit = 0.836 in the case where rj = and u = 1/q. The corresponding values of go are 
in the case of conservative mass transfer go = 3.77 for Pop. I giants (with parameters taken 
from Table 2) and go = 3.95 for Pop. II giants. In fact for all combinations (A, p) which 
describe the evolution of giants adequately we find go ^ 4. In the case of non- conservative 
mass transfer with 77 = and u = 1/q we find go = 1.89 for Pop. I giants and go = 1.97 for 
Pop. II giants, again with parameters taken from Table 2, and an upper limit of go ^ 2 for 
all acceptable combinations of (A,p). 

From the foregoing it is now clear that the near constancy of the mass transfer rate in the 
evolutions considered by WRS is not a generic property. Rather it is the result of a particular 
choice of initial parameters. In particular, if qi is close to (but larger than) g^it the mass 
transfer rate initially drops very strongly. On the other hand, the larger gf the more strongly 
the mass transfer rate increases beyond go. Because conservative mass transfer results in 
the largest possible values of gf we expect the mass transfer rate to increase substantially 
beyond go and that in particular if the initial mass of the primary is already much larger 
than that of the donor star. This is e.g. the case if the accretor is a black hole rather than 
a neutron star. To illustrate this we show in Fig. 5 the run of —M2 in a conservative mass 
transfer with Mi_i = 5M0, M2,i = 2M0 and Mc;i = O.26M0. It is clearly seen how the mass 
transfer rate grows because of the high final mass ratio of gf = 25.9. Nevertheless, the mass 
transfer rate stays reasonably constant if gi is not close to gcrit and the accretor is a neutron 
star, i.e. Mf i ^ 1.4M0. This is the case for the simulations shown by WRS. 

On the other hand, in the case that 1] = and u = 1/q the variation of the mass transfer 
rate is systematically smaller, first because for the same qi mass transfer is more stable 
than in the conservative case, and second because Mi remains constant, thus resulting in 
systematically lower values of gf. This is illustrated in Fig. 5 where we show the run of — M2 
for the same initial parameters as for the conservative evolutions shown in Figs. 3 and 4, 
but now with rj = and u = 1/q. As can be seen, the mass transfer varies much less than 
in the conservative case. 



9.2 Dependence on the initial core mass 



As can be seen from equation ( |32D , apart from the variation of the mass transfer rate during 
a mass transfer phase discussed in section |9.1|, the level of mass transfer is set by the initial 



core mass M^^i and depends on a large power, i.e. A — 1, of it. Using the parameters given 
in Table 2 which correspond to the results of WRS we have respectively —M2 oc Mc/ 
and —M2 oc Mc^^^ for Pop. I and Pop. II giants (see equations (|36|) and (^)). This strong 
dependence on Mc^ is a direct consequence of the strong dependence of the stellar luminosity 
on Mc,i (see equation (0))- Also the stellar radius shows a strong dependence on Mc,i as 
can be seen from equation (pIBI). Therefore, the initial core mass is the single most important 
parameter characterizing the evolution of a binary with mass transfer from a giant. 
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9.3 Dependence on 



We have already discussed the possible influence of Cc on the value of gcrit in section |9?1 
However, as Table 1 shows, the analytical solutions (0) - (0) depend on also via the 
exponent p2. Going from a solution with (^e = to one with (^e 7^ changes p2 by Ap2 = 
— Ccl-^^l) / P- Thus, with the parameters from Table 2 which correspond to the results of WRS 
we have Ap2 = — 7/5Ce for Pop. I, and Ap2 = — 4/3Ce for Pop. II chemical composition. 
Taking for an estimate Ce ~ —0.3 we find that for both cases, i.e. Pop. I and Pop. II, 
Ap2 ~ 0.4. Therefore, apart from increasing gcrit, going from = to Cc < leads to an 
increase of the mass transfer rate via the increased factor — Cr,2) and to a decrease by 
the factor (M2/M2,i)^^^. Sufficiently far from the stability limit, i.e. from gcrit, it is the latter 
factor which dominates. Because Ce ~ —0.3 this factor is never very small. For all practical 
purposes we find it to be between ~ 0.5 and ~ 0.7. Thus, qualitatively, the main effect of 
going from = to Co < is to increase gcrit and therefore to reduce the maximum allowed 
value for the initial secondary mass. 



9.4 Pop. I versus Pop. II systems 

Because for Pop. II giants Lq, Rq, A and p are significantly smaller than for Pop. I giants 
(see Table 2), the evolution of a Pop. II system, all other parameters being the same, is 
significantly slower than that of a Pop. I system if Mc^i ^ O.2M0. From equation (|32| ) it 
follows that the differential evolutionary speed scales as Mc,i^"~^' , i.e. using the parameters 
of Table 2 which correspond to the results of WRS, as Mc^"^. For an initial core mass of 
O.3M0 and O.4M0 the evolution of a Pop. II system is thus slower by respectively a factor 
of about 3 and 8. Furthermore, because Pop. II giants are also smaller than Pop. I giants if 
Mc ^ 0.2Mq the orbital period of a Pop. II binary is shorter too by a corresponding factor 
of about (Mc_i/0.2MQ)~^ than that of a comparable Pop. I system. 



10 CONCLUSIONS 

We have shown that the evolution of a binary with stable mass transfer from a giant donor 
can be described by a simple analytical solution (see section if the core mass-luminosity 
and core mass-radius relation are simple power laws of the mass of the degenerate helium 
core Mc. Furthermore, we have seen that the results of full stellar structure calculations of 
giants given by WRS and others can indeed be approximated with little loss of accuracy 
by simple power laws (section ^ and that, therefore, the analytical solutions derived here 
are relevant for the cases studied by Taam (1983) and WRS and similar cases. Moreover, 
the analytical solutions which we have derived are more general than the numerical ones 
of WRS in the sense that we allow for non-conservative mass transfer and a non-zero mass 
radius exponent Ce- 

A direct comparison of the results obtained from the analytical solution with those of 
numerical calculations by WRS (section |^) shows that both the analytical solution and the 
power law approximations are quite accurate (see Figs. 1-4). 

Using the analytical solution we have also made an a posteriori check of whether one 
of the key assumptions made in this model is fulfilled, namely that the donor star remains 
near thermal equilibrium despite being subject to mass loss (section |8.1| ), and have found 



that this is indeed the case. In section we put forward arguments, why, despite strong 
irradiation of the donor star from a compact accretor, the model is probably applicable to 
low-mass X-ray binaries. 
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Furthermore we have seen (section ^]1|) that the near constancy of the mass transfer rate 
over most of the mass transfer phase seen in the results by WRS is not a generic feature 
of this type of evolution but rather a consequence of a particular choice of parameters. In 
particular, we found considerable variation of the mass transfer rate if the final mass ratio q{ 
gets large. Therefore, particularly in conservative evolutions in which the accretor is massive, 
e.g. a black hole, the mass transfer rate increases quite substantially towards the end of the 
mass transfer phase (see Fig. 5). On the other hand, in evolutions in which the mass of the 
accretor remains constant or nearly so, qt gets much less extreme and as a consequence of 
that the mass transfer rate varies to a much lesser degree. 

From the analytical solution ( ^21) it is also directly seen that the level of mass transfer 
is largely set by M^^i, the core mass at the onset of mass transfer. 
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